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[57] ABSTRACT 

A method and apparatus for determining important histo- 
logical characteristics of tissue, including a determination of 
the tissue’s health. Electrical pulses are converted into 
meaningful numerical representations through the use of 
Fourier Transforms. These numerical representations are 
then used to determine important histological characteristics 
of tissue. This novel invention does not require rectification 
and thus provides for detailed information from the ultra- 
sonic scan. 

14 Claims, 5 Drawing Sheets 
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METHOD OF AND APPARATUS FOR 
HISTOLOGICAL HUMAN TISSUE 
CHARACTERIZATION USING 
ULTRASOUND 

5 

CROSS-REFERENCE TO RELATED 
APPLICATION 

This application is a continuing application of commonly 
owned patent application Ser. No. 08/592.833. filed Jan. 26, Q 
1996 and now U.S. Pat. No. 5,746,209. 

ORIGIN OF THE INVENTION 

The invention described herein was jointly made by 
employees of the United States Government and an 15 
employee of the University of Maryland Medical School. It 
may be used by and for the Government for governmental 
purposes without the payment of any royalties thereon or 
therefor. 
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BACKGROUND OF THE INVENTION 

1. Technical Field of the Invention 

This invention relates to the medical classification of 
human tissue as healthy or unhealthy. More specifically, the 2 s 
invention uses ultrasound for determining histological char- 
acteristics of tissue by converting the return energy pulses 
into numerical terms, thus facilitating a quantitative analysis 
for medical diagnosis. 

2. Discussion of the Related Art 30 

Ultrasound has routinely been used in the medical pro- 
fession to determine the shape, size and thickness of human 
tissue. Historically, this has been achieved through he detec- 
tion of specular reflections from macroscopic tissue inter- 
faces. Although this method allows for a general character- 35 
ization of the tissue, it does not reveal its underlying health 
and structure. In order to obtain such a detailed analysis, 
reflections from the microscopic interfaces (scatterers) con- 
tained within the tissue must also be detected. An improved 
invention to detect these scatterers has yet to be discovered, 40 

The prior art is able to determine general characteristics of 
tissue from its macroscopic interface reflections. This is 
commonly determined by transforming analog ultrasonic 
radio frequency data into a visual display. This technique 
requires the expertise of a skilled operator to interpret the 45 
display which does not include an accurate determination of 
the underlying pathology of the tissue. 

Other prior art systems use Fourier Transforms to detect 
microscopic reflections; however, these systems do not 50 
allow for direct measurement of reflectance through self- 
calibration, and therefore are subject to calibration prob- 
lems. These methods separate the components of an ultra- 
sonic pulse into corresponding frequencies, using the 
Fourier Transform. This requires a reference plate to nor- 55 
malize the Fourier energy from the reflected tissue and thus 
the accuracy and range of available data is limited. 

Similarly, the prior art taught by Sommer, Joynt. Carroll 
and Macovski ('‘Ultrasound characterization of abdominal 
tissues via digital analysis of backscattered wavefronts”, ^ 
Radiology , 141:811-7, 1981) uses a Fourier analysis to 
determine the mean spacing of scatterers in the liver and 
spleen. This method, however, does not provide for the 
determination of key variables that give ranges for specific 
tissues and their state of health. 65 

U.S. Pat. No. 5.417.215 of Evans et al. provides a method 
of interpreting the microscopic interfaces which requires 
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demodulating the return energy pulses by full-wave rectifi- 
cation to obtain the amplitude modulation of the pulses and 
computing the power spectrum by performing a Fast Fourier 
Transform on the rectified, digitized pulses. The requirement 
of rectification limits this method in providing information 
from the ultrasonic scan. It does not provide specific, well 
established physical parameters of the body tissue under 
investigation. It requires the comparison of Fourier energies 
or ratios, as well as. comparisons based on correlations. 

It is. therefore, an objective of the invention to provide an 
improved ultrasound tissue characterization system for his- 
tological tissue classification. 

It is a further objective of the invention to provide a 
system which does not require rectification, thereby allow- 
ing a more detailed analysis. 

SUMMARY OF THE INVENTION 

The present invention converts electrical pulses into 
meaningful numerical representations through the use of 
Fourier Transforms. This information can be used to deter- 
mine important histological characteristics of tissue, includ- 
ing a determination of the tissue’s health. This novel inven- 
tion does not require rectification and thus provides for a 
much more detailed information from the ultrasonic scan. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is an equipment diagram for data taking; 

FIG. 2 is an illustration of an ultrasonic wave passing 
through an iterated structure constructed from two acousti- 
cally different tissue types; 

FIG. 3 is a typical ultrasonic echo pattern received from 
tissue reflections; 

FIG. 4 is a plot of the Fourier Transform equation with 
arbitrary amplitude; and 

FIG. 5 is a plot of the Fourier Transform equation with 
discrete interpolation. 

DETAILED DESCRIPTION OF THE 
INVENTION 

Referring more particularly to the drawings, a well- 
damped 10 MHz transducer one-eighth inch (' «") diameter 
crystal 12 is connected to a wide band pulser-receiver 14. 
The pulser-receiver 14 is a sonic generator and receiver for 
generating, sending and receiving waves. The wide-band 
(non-rectified) output of the pulser-receiver 14 is connected 
to a digital oscilloscope 16. which has a plotter 18 connected 
to it to record the data output. The oscilloscope 16 provides 
live visual reading of the output. The data-point plotter 18 
gives a graphic reading of the output. The oscilloscope 16 
also has Fast Fourier Transform capability. 

The transducer 12 emits an ultrasonic pulse that travels 
through the body 20 into the subcutaneous tissue that is 
under examination. This tissue reflects a sequence of small 
amplitude echoes back to the transducer 12. The sequence of 
echoes from the tissue is received by the transducer 12, 
which converts the ultrasonic pulse sequence into a 
sequence of electrical pulses. The electrical pulse sequence 
is amplified by the wide-band amplifier in the pulser- 
receiver 14 and sent to the oscilloscope 16 for measurement. 
The oscilloscope 16, further processes the data, displays it 
on a visual monitor, and sends it to a plotter 18 for a 
hard-copy record. The oscilloscope 16 has the ability to 
perform Fast Fourier Transforms on the digitized pulses. 

Specifically, the transforming equations used are as fol- 
lows: 
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Gif) = 


1(a) 

hit) = 1 


1(b) 


where G(f) is the Fourier Transform of h(t). 

The following properties of this type of transform are 
used: 

F\k(t - *)>l = cxp{~i2nfto)F\h(t)] (shift operation) 2 (a) 

2(b) 


In general, any part of the transform may be used to 
express a relationship among the ultrasonic variables and, 
hence, give useful information about the significant tissue 
classifications. For purposes of illustration, the magnitude is 
used, although the real part, the imaginary part, the phase or 
any combination is just as instructive. 

Next, the Fourier transform of the incident ultrasonic 
pulse. 


F|A M (f)| (addition theorem) 

i»p=i I «=i 

F\ hi kt ) | = j g| j I ( similarity theorem) 
here t = yj - 1 


10 


15 


P - /’bexfl 


|cos2n/ 0 ib 


2(c) 
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Consider the echo train from an alternating stacked sequence 
of tissues. The following relationships are used: 

P 2 n— i = P o(T 1 7^2 1 ) 2 ^” — 1 'R 1 2 exp(~2n ,i,-2(/j-l) 2 / 2 )(odd mem- 
bers) 3(a) 

P 2 W =-Pofri 2 T 2 i) 2 ^ 1 Ri 2 exp(-2n| ,/,+ 2 / 2 ]Keven members) 3(b) 

In examining the echos, each received echo is subscripted 
with an index to indicate the temporal order of reception 


qj-n.-jco 

is evaluated by substitution into equation 1(a) 
F|/*o(0| = P,exp|-;r^ jcos(2*jS>f)</f 

which upon evaluation gives 


P O' rex p( _ m ~ fo I 

PIMrll = . / > 0 


( 8 ) 


(9) 


25 Using standard algebraic techniques and trigonometric sub- 
stitutions the magnitude is rewritten: 
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i -tflA/d* 

hjit - lb) = Pycxd — — — cos27r/ 0 (Ar/) 


(4) 


where ^ marks the reception of the first echo, and At^ is the 
time interval between the jth received reflection and the first 
echo. 

By allowing G/f^Flh^t-to) | and using (la), (2a) and (3a) 
the following relationship is obtained 


m m 

Gif) = Y Gj</)= Yj Flh)V-t))\=‘ifV r \k.W\ 

)=U. 5 >*2,4.6 

where K{f) = (^12 T^i y^ -1 Pi 2 ex p|"(ly + 4*1*1 - U-lfah) - 

j= I v 

* I I ll t 2 

y 1 Pi 2 «p| + (ilk] - *2/r//| — + — 


(5) 


- / 0 )) 2 } 

| Gif > 0, lb = 0)1 - X R n e\p - 

I 1 - 27cos(4/r ft 2 )exp< - 2H 2 ) + ( Texp< - 2 H 2 ) 9 
\ 1 - 27 2 cos(4^/r < )exp(-2//o) + (7 2 exp(-2A/o)) 2 


35 where f 0 is the center frequency of the ultrasonic transducer. 
Ho is otjli+a^ Hi is ajlj, H 2 is a^, T is (T l2 T 21 ), t 2 is 
lj/Cj, t, is lj/Cj+yc^ (Xl is the attenuation coefficient of the 
type 1 material, and a 2 is the attenuation coefficient of the 
type 2 material. 

An analysis of equation (10) shows the shape of IG(f>0. 
t=0)| is determined by two different undulations. The more 
rapid undulation occurs between much smaller frequency 
intervals, and is determined by the cosine function in the 
denominator of the radical. The cosine function, cos 47tft, 
causes this denominator to vary between the values 

[1-(T 2 exp(-2Ho))l 2 , and (1+<T 2 cxp(-2H 0 ))| 2 

both values are perfect squares. This occurs when cos 47tft^ 
varies between +1 and -1, for which 47rft =n7t or 


40 


45 


50 


*/,= 


(ID 


4-+-1 

\C 1 C 2 


It is assumed that the later reflections are small and equation where is the smaller frequency interval for the more 

(6) is not greatly affected by letting m approach infinity, the rapid undulation. The numerator inside the radical also 

equation is rewritten: varies between 


= 


Pi2 c *p|-2<*ifi - *4/r/~ 1 1 - iTnTu * X P -'hnk ~ i4^/~ 
|l -( 7 ’ l 2 7 ' 2 i) 2 expJ- 2 (c> l / l + (*2*2! ~ iAtn f ~ + ~)|] 


(7) 


65 


by substituting the sum for each series. 


| l-(T exp(-2H 2 ))] 2 , and [ 1+<T exp(-2H 2 )) j 2 
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but with a frequency of 


is used to obtain 


4 — 

<?2 


( 12 ) 


Z2 =Z t 


'V 

q 

€ - 1 ' 
,E + 1 , 

|exp( 2Hq i 

1+ \i 

q 

e - 1 ' 

U+ 1 

|exp<2// 0 1 


( 18 ) 


where Af^ is the frequency of the slower undulation. Solving 
equations (11) and (12) for 1 A and 1 2 respectively it is 
determined: 


10 


/1 

C\ 1 1 ^/« | 

(13a) 


“4A/.I A/2 | 


h 

C 2 

(13b) 


~4A/ 2 
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By measuring the frequency intervals, Af^ and Af 2 . one can 
determine an average thickness of a ground substance, 1 2 and 
the cluster of substantive cells, l x . 

The plot of equation ( 10) reveals a series of maximum and 
minimum values. One can trace the envelopes of the relative 
maxima and minima by noting the occurrence when cos 
471^=+! and -1, respectively 
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\OuppeAf > 0, tQ - 0)| = 


Ai r eXp(— r*{f - / 0 ) 2 ) 


(14) 


where Z x is the acoustic impedance of the cluster of fat cells, 
which is assumed to be the same as the acoustic impedance 
of the lipids stored in the fat cell clusters. These lipids 
possess a measurable acoustic impedance. Combined with 
the experimental determination of e, Z x and Z^ the 
acoustic impedance of the ground substance, can be calcu- 
lated. This analysis suggests a closer examination of the 
factors of the ground substance’s acoustic impedance. 

In determining the collagen and mass density |gms/cm| of 
tissues, the published values of the mass density of pure 
collagen varies from 1.16 to 1.33. However, collagen is only 
one of the many specialized molecular species found in 
tissues. Upon examination of the mass densities of the other 
constituents however it was discovered that a reasonable 
mass density model for ultrasonic applications is essentially 
a two-component mixture comprised of collagenous and 
other materials. Using the definition of mass density and 
further assuming that the collagen fiber diameters are 
smaller than the ultrasonic wavelength in the tissue it can be 
written 


R n ej.p,-2H l ) 


yjl- 27exp(-2// 2 )cosAtt fti + ( 7exp( - 2H 2 )? 


1 - 7^exp(- 2//o 1 


while 


\G lower {f > 0, t - O)) 


/? X2 exp(-2/y 1 ) 


A-expf-/r<^(/ - /o) 2 ) 


yj 1 - 27exp( - 2 Hi )cos4/r/f 2 + ( 7exp( - 2H 2 ))* 
1 + 7^exp(-2//o) 
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(19) 


C_ 

KX) 

Pi M 


C 

loo 


For a given frequency, f A (e.g. at a maximum for \G upp ^ i 
t=0)l) the ratio, 

_ \G uppe AfutQ=0)\ _ 1 + T 2 exp(-2// 0 ) 

€ \Gjower(fl , <0 = 0)| 1 - T^CXpi - 2 H 0 ) 


where C is the wet weight percent of collagen. p x is the 
density of the non-collagenous material in the tissue and p 2 
is the density of the collagenous material. 

The sound velocity in a mouse tail tendon fibers has been 
40 reported to be 1.733+0.056xl0 5 cm/sec (Goss and O’Brien. 
Journal of Acoustical Society of America. Vol. 65, pages 
./fj, 507-511; 1979). Tendon collagen is 30% wet weight. 

O’Brien (in “Proceedings Ultrasonic International 77”, IPC 
Science and Technology Press. Guildford, England, pages 
(16) 45 194-205; 1977) has empirically determined that the velocity 
of sound in tissue, c, depends upon the percent wet weight 
of collagen in the tissue’s composition. C, according to the 
following relation 


Solving for T 2 gives 

f£ - 1 \ 


50 




; xp( 2Hq ) 


(17) 
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where e is the ratio of maximum envelope to minimum 
envelope. Using data measurements one can obtain the 
values of e as a function of frequency and by curve fitting 
obtain values for T and H 0 . 

In determining the collagen content in the ground 
substance, equation (17) is solved for and 


c=( 1.58840.032 log*C)xl0 5 cm/sec 
The acoustic impedances are calculated to be 

(1.588+0.032 lo&Ox 10 5 , 

Z2 = — gins/ cur sec 

1 ~ loo + Too 

/>! 


( 20 ) 


(21) 


Zi -Z> 

R 12 = . R\2 — ~ R 21 

Zx+Zn 


The measurement and analysis technique gives Z 2 . With 
60 equation (21) it is possible to determine the weight percent 
of collagen in the ground substance in-vivo if the density is 
known of the collagen fibers. But the density of collagen 
fibers depends in part on the health and nutrition of the 
patient. This makes the value of Z 2 a potentially good 
65 indicator of tissue health with the larger values of Z 2 
indicating more and better grade collagen in the ground 
substance. By assuming a reasonable range of values for the 
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collagen fiber density the determination gives a range of 
values for the weight percent of collagen in tissues in-vivo. 
The linkage is the degradation of the mechanical properties 
of these tissues is directly involved with the failure of the 
tissue’s mechanical properties. That connection is the 
weight percent of collagen and the collagen’s density (the 
higher the density, the better its structural properties) in the 
ground substance. 

Using a curve fitting algorithm for the experimentally 
determined quantity 

| g ~ 1 ) 

Both T 2 and Ho can be determined by solving the equation 

/£-l\ (23) 

7 4 expt-2Ho) = (— y) 

where H is a 1 l 1 +a 2 l 2 . Since 1 A and 1 2 are determined 
according to equations 13(a). 13(b) and can be estimated 
from scattering and absorption of ultrasonic waves by fat 
cells, one can determine a 2 

Ho-ititi (24) 


which can be measured as a function of frequency. This 
information can be used to assess the quality of collagen 
fibers and the areal density of the fiber network (weave) in 
the ground substance, A correlation with invasive measure- 
ments and to a lesser extent to in vitro measurements will 
permit the best estimates of these quantities. 

In this manner, it is clear that more than a broad classi- 
fication of healthy or unhealthy can be realized for smaller 
regions of interest. Most significantly, a diagnosis for the 
development of pressure ulcers can be facilitated through the 
use of this invention. 

It is specifically intended that the present invention not be 
specifically limited to the embodiments and illustrations 
contained herein, but embrace all such modified forms 
thereof as come within the scope of the following claims. 

What is claimed is: 

1. A method for detecting histological properties of 
human tissue, composed of an iterated stacking sequence of 
cellular and collagen containing structural components, 
using data obtained from ultrasonic ray inputs, comprising 
the steps of: 

emitting a high-frequency ultrasonic pulse into human 
tissue; 

collecting a subsequent sequence of reflected unrectified 
echoes from said tissue; 

converting said unrectified echoes into a sequence of 
unrectified electrical pulses; 

amplifying said sequence of unrectified electrical pulses; 

performing a Fourier transform upon said sequence of 
unrectified electrical pulses; 

selecting any part of said transform; 

plotting a display of said transform in terms of said part; 
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evaluating said display to quantitatively determine a per- 
centage of collagen in said collagen containing struc- 
tural component, wherein said percentage can be used 
to assess the quality of said tissue. 

5 2. The method as in claim I wherein said selecting 

comprises using a magnitude value of said transform to 
determine histological tissue characteristics. 

3. The method as in claim 1 wherein said selecting 
comprises using a real value of said transform to determine 
histological tissue characteristics. 

4. The method as in claim 1 wherein said selecting 
comprises using an imaginary value of said transform to 
determine histological tissue characteristics. 

5. The method as in claim 1 wherein said selecting 
comprises using a phase value of said transform to deter- 
mine histological tissue characteristics. 

6. The method of claim 1 wherein said selecting step 
includes selecting any part of said transform to express a 
relationship among ultrasonic variables, and further com- 
prising the step of 

detecting pressure ulcers in muscle tissue based on said 
expressed relationship among ultrasonic variables. 

7. The method of claim 1 wherein said selecting step 
includes selecting any part of said transform to express a 
relationship among ultrasonic variables, and further com- 
prising the step of 

detecting pressure ulcers in adipose tissue based on said 
expressed relationship among ultrasonic variables. 

8. Apparatus for classifying histological characteristics of 
human tissue comprising: 

means for generating ultrasound pulses 
transducer means for applying said ultrasound pulses to 
said tissue and receiving returned ultrasound from said 
body 

means for converting said returned ultrasound pulses into 
a sequence of unrectified electrical pulses 
means for amplifying said pulses 
means for performing Fast Fourier Transforms on said 
amplified pulses 

means for determining said histological characteristics 
using any part of said transform. 

9. The apparatus of claim 8 wherein said means for 
determining histological charact eristics is adapted to deter- 
mine a magnitude value of said transform. 

10. The apparatus of claim 8 wherein said means for 
determining histological characteristics is adapted to deter- 
mine a real value of said transform. 

11. The apparatus of claim 8 wherein said means for 
determining histological characteristics is adapted to deter- 
mine an imaginary value of said transform 

12. The apparatus of claim 8 wherein said means for 
determining histological characteristics is adapted to deter- 
mine a phase value of said transform 

13. The apparatus of claim 8 wherein said means for 
determining is adapted to determine the histological char- 
acteristic of pressure ulcers in muscle tissue. 

14. The apparatus of claim 8 wherein said means for 
determining is adapted to determine the histological char- 
acteristic of pressure ulcers in adipose tissue. 
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